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Abstract. Two Component Advective Flow is the only complete solution 
that incorporates outcomes of actual theoretical solutions to explain spec¬ 
tral and timing properties of radiation emitted from the vicinity of black 
holes. It redefined the subject of black hole astrophysics by upgrading it 
from some sort of ’climatology’ and making it a precision science. Today 
any good spectral and temporal data could be fitted with TCAF with ease 
using as few as four parameters, totally un-heardof by the plethora of mod¬ 
els which are rat-racing to fit gross properties of data. TCAF addresses 
most of the issues of observations from both galactic and extra-galactic 
black holes while keeping the underlying framework (equations, assump¬ 
tions) unchanged. We discuss some of these points in this short review. 
As such, it concentrates on our group’s work to develop the subject till the 
present day. Most interestingly, these success of TCAF were accomplished 
without explicitly using any magnetic field. The magnetized disk solutions 
or simulated results in the literature till date are in the dark of how these 
issues could be addressed, indicating that magnetic fields in the accretion 
flow are either not implemented properly or may not be as important as 
they are thought and made out to be. Other models of hot accretion flow 
are either special cases of TCAF or are simply wrong and can explain some 
special features on special occasions. We claim that any observation that 
can be fitted by any of the multitude of models can be most certainly fitted 
with TCAF using far fewer number of parameters. 
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1. Introduction 


It is exactly twenty five years since it was realized (Chakrabarti, 1990ab, hereafter 
C90ab; 1996) that there is a fundamental change in topological properties of the 
accretion flow solution when viscosity of the disk is increased. In the language of 
a parameter prescription of Shakura & Sunyaev (1973), C90ab showed that when 
a < acritiE, A), a flow can have either one or two saddle-type sonic points. In the 
latter case, a standing or oscillatory shock can form. In the former case, the flow is 
hot (the exact version of the all the present ‘hot flow models’) and passes through one 
sonic point very much like a Bondi flow (1952). However, when a > acrit, there is 
only one sonic point, and disk is sub-sonic and Keplerian beyond this point. A quote 
from 1993 conference proceedings is in order: “These findings are very significant 
as they propose a unifying view of accretion disks. This incorporates two extreme 
disk models into a single framework ..” (Chakrabarti, 1994). This dichotomy of flow 
topology is the genesis of the two component advective flow (TCAF) solution which 
combines the latter type solution (a > acnt) on the equatorial plane and the former 
type solution (a < acnt) away from the plane (Chakrabarti, 1995a; Chakrabarti & 
Titarchuk, 1995, hereafter CT95). 

Prior to RXTE, an understanding about X-ray spectral states of stellar mass black 
holes has really not taken a concrete shape. More studies were available on Active 
Galactic Nuclei. Not surprisingly, initial models of TCAF were dedicated to super- 
massive black holes (Chakrabarti, 1994, 1995ab, Molteni, Sponholz and Chakrabarti, 
1996). CT95 was written when RXTE was still being fabricated and tested anticipat¬ 
ing that precision science on stellar mass black holes could be possible with RXTE. 
While TCAF was being sporadically used for fitting some RXTE data (Chakrabarti et 
al., 2005; 2008; Chakrabarti, Dutta & Pal, 2008), it was not until numerical simula¬ 
tions of Giri (& Chakrabarti (2013) where it was established that TCAF was a stable 
configuration, XSPEC modeling with TCAF was tried out. A few preliminary papers 
have been recently written to establish that CT95 model and its minor improved ver¬ 
sions based on rigorous theoretical solution is not only capable of explaining all the 
issues related to the spectral properties and extract the most relevant physical para¬ 
meters, it also explains majority of the timing properties, such as evolution of these 
parameters. It is a promising tool for future missions, for black holes of all masses. 
The effects of spin will be incorporated in TCAF in due course. 


2. TCAF 


A cartoon picture of TCAF was presented in the first RETCO conference (Chakra¬ 
barti, 2013). Basic properties of TCAF are also discussed there. Briefly, the compon- 
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ent with higher viscous shear stress becomes Keplerian and if accretion rate is signi¬ 
ficant, it produces standard disk and emits soft photons by efficient cooling. These 
are intercepted by the advective component of lower viscosity, especially the region 
where standing or oscillating shocks (CENBOL) are produced. When shocks are ab¬ 
sent, the flow will still be denser due to the centrifugal barrier (Chakrabarti, 1997), 
so even the hot accretion flows (transonic flows passing from outer sonic points) can 
be useful. The intercepted photons are inverse Comptonized to produce the ‘power- 
law’ component (C95). When the CENBOL oscillates, either due to a resonance 
between the cooing time scale and the infall timescale (Molteni, Sponholz & Chakra¬ 
barti, 1996; Giri, Garain & Chakrabarti, 2015; Chakrabarti et al. 2015), or because 
Rankine-Hugoniot condition is not fulfilled (Ryu, Chakrabarti & Molteni, 1997) one 
observes quasi-periodic oscillations. These QPOs are sharper (high ’Q’, type C or 
B; Casella et al. 2005). In a hot flow with the centrifugal barrier but no-shock, the 
oscillation is not sharp, but some low Q-factor (Type A) QPOs are possible. Recently, 
spectra of some black holes using TCAE were fitted and disk/halo accretion rates, 
shock location and shock compression ratio (by these we mean typical length scales 
of the centrifugal barrier and its strength even when shocks are absent) were extrac¬ 
ted. This fitting routine is being improved (Debnath et al. 2012, 2013, 2014, 2015ab, 
Mondal et al. 2014, 2015) to take care of broader spectral properties. Pitted paramet¬ 
ers from spectra also enable one to explain timing properties, such as QPO frequency, 
a feat only achieved by TCAE Chakrabarti et al. (2015) shows that LPQPOs may 
form only when Tr, the ratio of the cooling time scale (C to infall time scale f,) is 
around unity (say, when 0.5 < < 1.5), 


T, = - = 3.5 X 10^ 


^+Ar 


-d 


R2' 


where, Ar = {Mh + Md)IMd is the ratio of the total accretion rate to disk accretion rate, 
/o is the fraction of soft photons intercepted by the CENBOL, E is the average energy 
enhancement per soft photon after going through the CENBOL and R is the ratio by 
which inflow velocity of the advective component drops at the CENBOL boundary. 
This is a mass independent expression and should be valid for even super-massive 
black holes. In an outburst source, in hard and hard-intermediate states, satisfies 
resonance condition in general, and LPQPOs caused by resonance are seen in these 
states. In soft intermediate states, the flow is far from resonance. Its LPQPOs are due 
to non-compliance of Rankine-Hugoniot shock conditions even when the flow has 
two saddle type sonic points. In soft states, usually there is one (inner) saddle point 
and thus none of these QPOs occur. Only high frequency oscillations of the inner 
sonic point cause high frequency QPOs. 


3. TCAF vs. other disk models 


Since every feature of TCAE arises out of ‘some’ theoretical solution or numerical 
simulation, it is difficult to compare with other models which are ad hoc, empirical 
or phenomenological. They require a large number of free parameters (usually 6-10 
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or more) for any fit. Even normalization constant which is supposed to depend on 
distance, inclination and mass of the object is made variable! There are piece-wise 
concrete solutions such as: Bondi flow (Bondi, 1952), thin disk (Shakura-Sunyaev, 
1973), thick disk (Paczynski-Wiita, 1980; Chakrabarti 1985 and references therein) 
and generalized Advective disks with and without shocks (Chakrabarti, 1996). TCAF 
is a combination of all these, in general, and much more. In this scenario, hot solutions 
such as ADAF/RIAF etc. would be shock free component of Chakrabarti (1989 or 
1996) solution when (invalid) assumption of self-similarity is imposed on. Indeed 
workers get better solutions of ADAF following our procedure: “Chakrabarti and his 
collaborators (see Chakrabarti 1996for references) introduced a very clever mathe¬ 
matical trick... In this way, the most difficult part of the problem aA§ finding the 
eigenvalue dA§ is trivially solved.. " (Igumenshchev et al. 1998) or "Chakrabarti 
and his collaborators introduced a very clever procedure (e.g., Chakrabarti 1996a). 
The difficulty of finding the eigenvalues was simply avoided" (Fu, Gu & Yuan, 1999). 
Evidently, The so-called ADAF also has shocks (Beker et al. 2011) making these 
‘hot flows’ special cases of C90ab and Chakrabarti (1996) solutions. Apparently, any 
observation where the radiation flux does not follow the standard disk (SS73) value is 
assumed to have ‘ADAF’. Even a ‘slim-disk’ is also ‘ADAF’! 

Since CENBOF in TCAF behaves as a boundary layer (Chakrabarti et al. 1996), 
jets in black hole candidates are originated at CENBOF. The ratio between jet and 
disk rates (Chakrabarti, 1999; Das &. Chakrabarti, 2000; Das, Chattopadhyay, Nandi 
& Chakrabarti, 2001; Singh (& Chakrabarti, 2012; Kumar et al. 2013) are computed 
easily. TCAF thus includes ‘Jet dominated’ solutions. CENBOF concept is used 
in Jet-ADAF solution (e.g.. Yuan et al. 2002ab, Markoff 2003). M87 jet is already 
known to be ejected from a few tens of Schwarzschild radii (Junor, Biretta & Fivio, 
1999) justifying our conjecture. As Garain et al. (2012) showed, the jets are quenched 
when the base of the jet is cooled down (e.g., in a soft state). This type of prediction 
is out of reach of any other model. TCAF can incorporate jet emission as it is an 
extension of the CENBOF itself. 

TCAF does not yet explicitly include magnetic fields, as it is not clear if large 
scale magnetic fields could be sustained around black holes or if the field plays any 
major role. Highly magnetized plasma is notoriously unstable. Indeed, to our know¬ 
ledge, no single GRMHD code, includes the plasma instability or ever found jet form¬ 
ation with random poloidal fields. The field configuration has always been carefully 
choreographed. The simulations do not run beyond a couple of dynamical timescales. 
So the results are not trustworthy. It is doubtful if a magnetized disk can exhibit 
precisely the same type of variabilities as seen in, e.g., in GRS 1915H-105 after long 
time intervals. However, in 6 and [i classes (Belloni et al., 2000), toroidal magnetic 
field may play a major role in sudden disappearance of CENBOF, creating soft states 
and radio flares (Nandi et al., 2001). The plasma instability which causes sudden 
disappearance of CENBOF may also be responsible for production of super-fast jet 
components which are then collimated by predominantly toroidal flux ejected from 
disks due to buoyancy. Indeed, no fitting model in the literature uses any of the MHD 
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solutions which exist in the last four decades. On the contrary, as soon as stability 
of TCAF configuration was established, it was immediately implemented in XSPEC 
and is routinely used to fit spectral and timing data. Extraordinary insightful flow 
parameter values are extracted with surgical precision. Very good understanding of 
evolution of the accretion rates of the two components during an outburst has emerged 
from these fits. Success of TCAE has not gone unnoticed. 


4. Predictability of TCAF 


Today there are more models than the number of objects to fit and new names of 
accretion models are coined every day! While each such model was advanced with a 
specific goal to fit certain aspects of observations, the theory of transonic astrophysical 
flows (C90b) had no specific goal in mind. It was created to seek the most general flow 
solution in a black hole geometry. As such, all other disk models, if correctly solved, 
must be a special case. TCAE is a combination of two most important topologies of 
these transonic flows. 

So far, good fit of a given data is all that is required to be a good model - no link¬ 
age between successive observations of the same object or one object and another. A 
weak linkage was sought (in diskbb + Power-law model, for instance) through vari¬ 
ation of normalization constant which apparently gives which often turn out to be 
less than the marginally stable radius, or even less than the horizon! This is absurd. It 
arises out of misconception that normalization of an entire spectra should have some¬ 
thing to do with just the inner edge of the Keplerian flow which affects only a limited 
energy range. To this effect, TCAE has totally done away with normalization constant 
variation. Successful fits are made for a given object with very limited variation of 
the normalization constant at all (Molla et ak, 2015, in preparation; Jana et al. 2015, 
submitted). 

TCAE directly addresses the cause of LEEQPOs and predicts quite naturally how 
it should evolve. In outbursting sources, as the Keplerian accretion rate is increased 
due to increase of viscosity, the shock moves inward (Molteni, Sponholz (& Chakra- 
barti, 1996; Garain, Ghosh & Chakrabarti, 2014; Mondal et al. 2015) and QPO fre¬ 
quency rises. Once resonance condition sets in, it is locked over a range of flow 
parameters (Chakrabarti et al. 2015). Eitting of spectra during an outburst by TCAE 
shows this relation between shock location and QPO frequency. 

Then there are issues about color-color diagrams which TCAE squarely addressed. 
Historically, these diagrams were drawn for neutron stars and workers confidently 
talked about Z-sources, Atoll sources, banana sources etc. (e.g., Schultz et al. 1989). 
Blindly transporting these concepts to black hole physics would not help since stellar 
black holes could have masses from 3 to ~ 3QMq, unlike the case of neutron stars 
where the masses are always within a narrow range. With mass, the meaning of soft 
and hard photons changes and the shape of color-color diagram would be different for 
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different mass of the black hole or the accretion rate. (The very reason why hardness- 
intensity diagram or the so-called ‘q-diagram’ for outbursting black holes rarely look 
like ‘q’! See, Chakrabarti, 2008). TCAF, by having the Compton cloud as a part of 
the disk itself, demands that the mass of the black hole directly decides the sizes of the 
soft and hard photon sources. As a result, the so-called Comptonizing efficiency (CE), 
i.e., hardness ratios with soft and hard photon range dynamically determined (Pal & 
Chakrabarti, 2013), is a notion independent of mass. No wonder then that the CEs 
of GRS 1915H-105 and IGR 17091-3624 are found to be identical for a given variab¬ 
ility class, establishing a much deeper aspect of accretion physics that the geometry 
variation in these two objects are of similar nature. In fact, it also predicts that the 
appearances of variability classes would be identical in both the classes. Time scales 
being grossly different (heartbeats of IGR in the so called p class of IGR is faster by a 
factor of three as compared to that of GRS), it is unlikely that the masses are similar 
(Pal & Chakrabarti, 2015; this volume). TCAE is also an excellent tool to obtain mass 
of black holes (this volume). 

In binary systems, where properties of accretion flows, including spectral class 
and timing properties, are dependent on boundary values, definite prediction is diffi¬ 
cult unless the evolution of the boundary conditions, i.e., the pattern of mass supply 
from the donour star is known well known ahead of time. In the context of galax¬ 
ies, long term prediction may be possible since the variability time scale itself is very 
large. 


5. Whither TCAF? 


TCAE faces two common questions: (1) ‘Why the shock was introduced at all? Was 
it necessary?’ The question is meaningless, since shocks were never ‘introduced’. It 
is a solution of any valid set of equations of test flow around a black hole for a large 
range of boundary conditions. It is the direct consequence of having a horizon which 
is responsible for the innermost sonic point. Eor a neutron star, the innermost sonic 
point is inside the object and thus the shock forms near the boundary (the well known 
boundary layer). So this shock has as much right to be there as the boundary layer of 
any star. Indeed, it is shown recently that the highest Ucrit for which a shock can still 
form is ~ 0.3 (Nagarkoti & Chakrabarti, 2015a; this volume) though parameter space 
shrinks significantly. Since viscosities obtained by numerical simulations are around 
0.01 - 0.05 (Parkin, 2014; Sano et al. 2004), the fact that shocks should exist in many 
cases is a foregone conclusion. However, outside of this parameter range with two 
saddle points, oscillating shocks would still be present causing QPOs. In extreme 
cases, a low viscosity flow would always have an excess centrifugal pressure which 
will slow down matter. Observationally, presence of Type C QPOs (Caselli et al. 
2005) are evidences that the Comptonizing region has a sharp boundary. I conjecture 
that at the outskirts of resonance condition. Type B can form. And when the shock is 
not involved, and only centrifugal barrier is involved. Type A QPOs with very poor 
rms and low Q values cannot be ruled out. 
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(2) ‘how does one supply sub-Keplerian flows?’ First of all, to maintain a Kep- 
lerian distribution in a disk, a significant viscosity is required (a > acrit) which, ac¬ 
cording to recent simulations (Sano et al. 2004; Parkin, 2014) is hard to come by. It 
is difficult if not impossible to remove 99.82% of specific angular momentum from 
a Keplerian injected matter (say, from ~ 1000 at outer edge to ~ 1.8, the marginally 
stable value) within an infall time scale (of even a couple of hours in some of the 
MAXI sources). It is much easier to redistribute low angular momentum flow into 
a Keplerian disk by moderate viscosity. It is not surprising that objects are found to 
stay in soft states for very short periods as compared to their duration in hard states. 
So this low angular momentum component must exist. Most interestingly, the same 
inquisitive workers use imaginary hot flows of unquantified property and shape very 
often when a standard Shakura-Sunyaev disk is not found adequate. In a more seri¬ 
ous note, our low-angular momentum component is natural for any high mass X-ray 
binary system (Smith et al. 2001, 2002), and most certainly for Active galaxies where 
matter is supplied from stellar winds (Chakrabarti, 1994). Again, not surprisingly, 
for M87, when TCAF was used to fit the optical line emission data, the mass was 
found to be ~ 4 X 10^ Mq (Chakrabarti, 1995b), closer to the currently accepted value 
(Walsh et al. 2013) of 6.6 x IO^Mq than 2 x IO^Mq, assuming SS73 (Harms et al., 
1994). Since precise nature of black hole spectra cannot be left to poorly understood 
properties of corona, evaporated disk etc. the sub-Keplerian matter has to be supplied 
partially even for low-mass X-ray binaries, even if it means that matter has to loss 
most of its angular momentum on the donor surface itself through coronal fields. 
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